不同围压作用下非均匀岩石水压致裂过程的数值模拟

从岩石细观非均匀性的特点出发，提出一个描述非均匀材料渗流和破裂相互作用的数值模型。在这个数值模型中，单元的力学、水力学性质根据统计分布而变化，以体现材料的随机不均质性，材料在开裂破坏过程中流体压力传递通过单元渗流，损伤耦合迭代来实现。算例表明，该模型能较好地模拟出岩石类材料在水力压裂作用下，微结构非均匀分布和不同围压比对破裂模式、失稳压力的影响，非均匀性导致试件的开裂压力、失稳压力明显不同，裂纹扩展路径不规则发展，模拟结果和实验结果较为一致。     

心电激励下的人体左心室力学响应分析

人体心脏通过电兴奋引起的心肌收缩实现泵血功能 ,而心肌的力学特性高度依赖于肌纤维结构。本文根据肌纤维旋向和复合材料理论以及电生理心脏模型建立了左心室的有限元机械模型 ,仿真研究了左心室在心电兴奋力作用下的力学响应。结果表明左心室的收缩过程是十分复杂的 ,包括轴向和径向的收缩及绕长轴的不同程度的旋转扭曲。从总体趋势上看 ,心尖处变形最为严重。此外 ,仿真结果还表明 ,心壁应力分布不均匀 ,从内壁到外壁有所减少 ;在心尖和心底部应力较大 ,其中内壁心尖处应力最大。这些结果说明了心脏的力学特性与心肌纤维结构、左心室的几何形状以及电兴奋刺激密切相关。作者将本文的仿真结果与医学图像及其它模型进行了对比分析 ,结果表明了该模型用于分析心肌力学功能特性的可行性     

FV/MC混合算法求解轴对称钝体后湍流流场

介绍一种有限容积／Monte Carlo结合求解湍流流场的相容的混合算法．有限容积法求解Reynolds平均的动量方程和能量方程，Monte Carlo方法求解模化的脉动速度—频率—标量联合的PDF方程．将该算法发展到无结构网格，探讨了在无结构网格中实现两种方法的耦合，包括颗粒定位，颗粒场和平均场之间数据交换等问题．并以二维轴对称钝体后湍流流场作为算例，比较了计算结果与实验结果．     

Numerical simulation of Mach reflections

In this paper, the “FLIC” difference method with triangular mesh is adopted to numerically simulate the regular and Mach reflections that occur when a shock wave pass around a wedge. The compuational result is compared with the shock tube experimental results of G. Ben-Dor and I. I. Glass. The comparison shows that the position, shape of shock wave and height of Mach stem all show a good agreement. Consequently, the “FLIC” difference method with triangular mesh is quite satisfactory in numerical simulation of the regular and Mach reflections.     

Analysis of temperature and stress field in Al alloy’s twin wire welding

By analyzing the shape of twin wire welding’s arcs and the track of droplets’ transition, the phenomenon that the twin wire welding’s fore arc and rear arc all deflect to the middle of the two arcs is found. Based on this the double ellipsoid heat source model is amended, and a heat source model which can apply to calculate the twin wire welding’s temperature field is put forward. This model is testified by actual experiment of temperature sampling. By comparing the temperature field of twin wire welding and single wire welding, the results show that twin wire welding has slender weld pool the end part of which is ellipsoid, and its HAZ is narrower than that of single wire welding. So, twin wire welding can not only reduce the Al alloy generating hot crack, but can also weaken the “overaging” softened phenomenon of heat treated strengthening Al alloy. In the end, the evolving rules of 2219 Al alloy’s longitudinal and transverse stress when welded with twin wire welding are analyzed.     

平面叶栅跨音流场气动特性的数值研究

通过对模型方程的分析，给出了一种新的隐格式构造思想。将它运用到关通量分裂格式中，可得到无近似因子分解、无矩阵运算的高效二阶精度隐式矢通量分裂差分格式，并用来直接求解时间平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组。数值计算标明：该方法具有精度高、稳定性好、计算量少、收敛快等优点，在平面叶栅跨音流场的计算中，较好地捕获了激波，与实验比较，结果令人满意。     

逆向环形分布射流迴流特性的数值计算

数值计算表明,逆向多环形分布射流可造成颗粒能直接进入,并能进行强烈的动量、能量交換的理想的迴流区结构。多环形结构可在大范围內造成所需尺寸的迴流区。逆向射流和主流速度比、环直径、射流直径和孔数、环的结构形式及离燃烧室进口的距离是影响迴流特性的主要参数。     

挤出平缝口模优化设计

挤出平缝口模通常用于加工膜和片材，对产品厚度的一致性有很高要求。本文给出了结合聚合物成型模拟技术、设计灵敏度分析和数值优化技术的口模形状优化设计方法。以最小压力降为目标函数，口模出口处任意点的速率与出口已知速率相对误差的平方和小于容许误差为约束函数，口模形状参数为优化设计变量，给出了目标函数的表达式，推导了日标函数、约束函数对优化设计变量的灵敏度公式。利用灵敏度分析和基于梯度的优化算法即序列二次规划算法(SQP法)求解最优设计参数。通过算例表明，应用该法进行口模优化设计可以减小压力降和口模出口速率变化率。     

粘塑性损伤模型模拟准超塑性单轴拉伸行为

发展了Chaboche粘塑性本构模型的大变形隐式算法，用损伤(DM)和无损伤(NDM)模型模拟准超塑性单轴拉伸。发现变形过程可分为三个阶段：均匀变形、颈缩发展、断裂破坏阶段。DM可准确模拟前两个阶段变形，NDM只能较好地模拟均匀变形阶段，表明DM可以较精确地描述稳定发展的动态过程。由于有限元方法只能描述连续介质，因此对于断裂破坏阶段，NDM模拟载荷大于试验结果，DM的载荷小于试验结果，这是由高应变速率敏感性造成。DM能够描述试验中出现地多处颈缩现象，局部应变速率分布随时间演化反映了颈缩发展程度。严重颈缩部位的距离代表着超塑性变形能力，距离越大，抗颈缩能力越好。